Erratum: Initial fields and instability in the classical model of the heavy-ion collision 
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We correct a mistake in the analytical expression for the energy density given in Phys. Rev. C 76, 
021902 (2007) _arXiv:0704.3625 [hep-ph]]. The expression should be multiplied by 16. One question 
then arises; how could it be possible to explain this difference between the analytical and numerical 
results in the same model if both are correct? We find a subtle problem in the treatment of the 
randomness of the color source along the longitudinal direction and the treatment of the longitudinal 
extent of the color source. 
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The initial energy density given in Ref. 1] should be 
multiplied by 16, that is, Eq. (6) should be corrected as 



2tr{Bi;)') = \N,{Nl-l)a\, 



and, accordingly, Eq. (8) should be corrected as 
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Also, not only the r^-order terms but the r^-order terms 
missed the same factor 16. Equation (11), therefore, 
should be 



and Eq. (13) should be 
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As a result, the estimate for the initial energy density 
in Eq. (14) (which had a typo; tt in the square brack- 
ets should be ^), its resummed form in Eq. (15), and 
the compact formula in Eq. (16) in the continuum limit 
should be corrected respectively as 
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obtained from Af^''^ derived as Eq. (75) in Ref. [31]), while 
Eqs. (16) in Ref. 0] seems to be based on Eq. (2.23) in 
Ref. [J|. These are independent calculations with differ- 
ent cut-off prescriptions, but of course, the coefficients in 
front of the logarithmic singularity should coincide with 
each other. 

In a quantitative sense, in fact, the analyses on the 
Glasma instability should be affected by the missing 
factor 16 through Eqs. (22) and (23) in Ref. [l|. Our 
discussion on instability is not beyond the qualitative 
level, however, and the essential idea for the possible 
instability mechanism still works. Also, in the first part 
of Ref. 1], no modification is required in the essential 
ideas; the energy density at r = has a logarithmic 
divergence and the expansion in terms of r is singular 
around t = 0. The resummed form in the logarithmic 
ansatz is a natural choice to remove the singularity at 
finite r. 

One question arises immediately, however. The initial 
energy density has been evaluated in the same model and 
the same cut-off prescription in different two methods. 
If both methods equally work well, two results must be 
identical apart from a small discrepancy originating from 
the lattice and continuum formulations. Why can one 
differ from the other? 

This problem is so interesting on its own that it may 
deserve another publication but we shall briefly see 
where the difference stems from. In short, the important 
point is that the McLerran-Venugopalan (MV) model im- 
plemented in the numerical calculation is not faithful 
to the analytical formulation. The subtlety originates 
from how to define the ill-defined expression of Eq. (4) 
in Ref. We should introduce some regularization to 
write it in a form of 
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and 
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We remark that the first term of above Eq. ([5]) has the 
same overall factor as given in Eqs. (14) and (16) in 
Ref. [1] . Our calculations make use of the technique 
developed in Ref. Q (the necessary correlation function is 



Here we defined the regularized color source as 
lim p^{x±,x^) = p{x±) S{x^) . 
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Another delta function in the longitudinal direction ap- 
pears in the correlation function in terms of the sources, 
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where we replaced the delta function by the regularized 
one in the longitudinal direction such that, 



To evaluate the left-hand side, we have to perform the 
following Gaussian integral. 



lim 5c(x ) ~ S(x ) . 
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The question we are addressing here is whether we are 
allowed to adopt the following simplification; 
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to prevent the delta function from appearing at all. The 
numerical calculations commonly make use of Eq. for 
practical reasons with hope that the final answer would 
not depend on this replacement (see, e.g. Eq. (36) in 
Ref. 0, Eq. (4) in Ref. ;6], and so on). 

Using the notations we introduced above, we can re- 
iterate the question more rigorously. That is, we shall 
check, 



hm\un{0[V,]). = hm lim(0[y,]) 
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The left-hand side corresponds to the numerical imple- 
mentation and the right-hand side to the analytical for- 
mulation of the MV model. 

We treat more general cases in a separate literature [7] 
and will limit the current discussion only to the tadpole 
calculation, namely, 0[V] = . This simplest choice is, 
as we will see, enough to exemplify two limits in Eq. (jl3p 
are noncommutative. 

We already know the analytical answer for the right- 
hand side of Eq. (O. That is given by [1, H 



lim lim(vj^)^ — exp 
where fi^ — J dx^ fj,'^ {x^ ) and 
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where i°'s are the SU(A^c) algebra. It is hard to do this 
integration for arbitrary SU(A^c) group, while the SU(2) 
case is immediately doable. 

In the case of SU(2) (i.e. Nc — 2), the above Gaussian 
integral leads to 
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which is obviously different from the right-hand side with 
Nr ~ 2 substituted, that is, 
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from Eq. (|T4l) . 

Interestingly, though Eqs. p7|) and p8)) look quite dif- 
ferent, the curvature near 77 is the same. As we 
report in Ref. 0, the discrepancy between the left-hand 
and right-hand sides in Eq. (|13p becomes significant when 
the singlet component survives unsuppressed by the sys- 
tem size. In fact, the expectation value of the gauge 
fields, {VediVj ■ VediVj')(, leads to different energy den- 
sities by a factor. 
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